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摘要 :物种 地 理 分 布 主要 取决 于 它 对 气候 、 地 形 等 环境 因子 的 适应 性 。 基 于 22 个 环境 因子 和 阔 叶 红 松林 的 4 类 主要 建 群 树 
种 一 一 红 松 、 紫 机 .水 曲 柳 和 蒙古 栎 的 地 理 分 布 数据 ,采用 最 大 箭 模型 模拟 了 了 疮 叶 红 松林 的 潜在 分 布 区 域 , 并 分 析 决 定 阔 叶 红 松 
林地 理 分 布 的 主要 气候 和 地 形 因子 ,最 后 利用 政府 间 气 候 变化 专门 委员 会 (了 PCC ) 发 布 的 3 种 排放 场景 (SRES-A2、SRES-A1B、 
SRES-B1) 下 2020 ,2050 .2080 年 的 气候 数据 预测 阔 叶 红 松林 的 未 来 潜在 分 布 区 。 结 果 表 明 : 各 树种 的 受 试 者 工作 特征 曲线 下 面 
积 (AUC 值 ) 都 大 于 0.8 ,说 明 模型 有 很 好 的 预测 能 力 ; 影 响 阔 叶 红 松林 分 布 的 主导 环境 因子 是 年 降雨 量 . 季 节 人 性 降雨 量 海拔 、 
年 平均 温度 .最 湿 季 度 的 平均 温度 。 在 基准 气候 条 件 下 , 阔 叶 红 松林 的 高 度 适宜 分 布 区 主要 分 布 在 长 白山 和 小 兴安 岭 地 区 , 占 
研究 区 总 面积 的 11.69% , 低 度 适宜 区 面积 不 适宜 区 面积 分 别 占 研究 区 总 面积 的 23% 和 65.31% 。 模 型 预测 结果 显示 ,未 来 在 
A2 AI1B 和 B1 气候 情景 下 , 阔 叶 红 松林 高 度 适 宜 区 的 南 界 与 北 界 都 向 北 移动 ,其 面积 有 缩减 的 趋势 ,而 低 度 适 宜 区 的 面积 有 增 
加 的 趋势 。 

关键 词 : 阔 叶 红 松林 ;气候 变化 ;最 大 灶 模 型 ;潜在 地 理 分 布 ;气候 变化 情景 
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Abstract: Species-environment relationships have always been a central issue in ecology and biogeography. The adaptation 
of species to a changing natural environment, including topographical and climatic factors, determines its geographic 
distribution. Accumulating evidence has revealed that Northeast China has experienced the greatest temperature increase 
since the 1950s, which in turn affected species distribution in this region. Broadleaved Korea pine (Pinus koraiensis) mixed 
forest ( BKF) is the native forest type in the southern part of Northeast China. The distribution area of this species has 
shrunk substantially due to historical overexploitation. It is important for forest managers to be able to predict the potential 
geographic distribution of BKF based on species-environment relationships. Although the distributions of individual species 
in BKF such as 已 koraiensis and Quercus mongolica have been reported previously, few studies have focused on the 
potential geographic distribution of BKF. In this study, four dominant tree species -P. koraiensis, Q. mongolica, Tilia 


amurensis and Fraxinus mandshurica, which together account for more than 80% of the growing stock in primary BKF, were 
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selected to represent this forest type. Nineteen climatic and three topographic variables in Northeast China that are 
considered to be most likely influences on the geographic distribution of tree species were selected as environmental factors. 
To identify the major climatic and topographic factors controlling BKF distribution and simulate the potential geographic 
distribution of BKF under current climatic condition, the geographic distribution records of the dominant tree species, 
together with the environmental factors, were used in the Maxent model. The future geographic distributions of BKF were 
consequently predicted for the 2020s, 2050s, and 2080s, based on three kinds of climate change scenarios ( SRES-A2, 
SRES-A1B, SRES-B1) published by the Intergovernmental Panel on Climate Change (IPCC). For the four dominant tree 
species, the simulations showed that the area under the curve indexes ( AUC) were 0.925, 0.890, 0.859, and 0.847, 
respectively. All these values exceeded 0.8, which indicates that the models had a good predictive performance. The major 
environmental factors affecting the distribution of BKF included annual precipitation, precipitation seasonality, elevation, 
annual mean temperature and mean temperature of wettest quarter. For the entire region, 11.69% of the total area was 
identified to be of high suitability for BKF distribution, 23% was of low suitability, and 65.31% of the area was unsuitable. 
Under the A2, AlB, and B1 scenarios, the model predicted that both the southern and northern boundary of the high 
suitability area for BKF will shift northward. Overall, the high suitability area in this region was predicted to decrease, with 
the extent of the decrease depending on the severity of climate change. For example, under the A2 and AlB scenarios, the 
high suitability area in Northeast China will be less than 1% of the region by the 2080s. Overall, the results indicate that if 
no effective measures are taken to mitigate climate change, there is a great possibility that BKF will disappear from 


Northeast China. 


Key Words: broadleaved Korea pine mixed forest; climate change; Maxent model; potential geographic distribution; 


climate change scenarios 


分 析 物 种 -环境 之 间 的 关系 已 经 成 为 生态 学 .生物 地 理学 中 的 一 个 焦点 问题 1 ,物种 在 各 种 各 样 的 生物 因 
子 和 非 生物 因子 的 共同 作用 下 ,都 有 其 暂时 的 地 理 分 布 范围 。 物 种 地 理 分 布 主要 取决 于 它 对 气候 地 形 等 环 
境 因子 的 适应 性 ,特别 是 森林 群落 优势 物种 的 地 理 分 布 与 气候 有 着 密切 的 关系 ,其 地 理 分 布 受到 气候 变化 的 
深刻 影响 后 。 而 很 多 大 气 环流 模型 预测 研究 表明 ,到 21 世纪 末 ,地表 平 均 温 度 将 升 高 1.1 一 6.4%5 ,因此 ,全 
球 气候 变 暧 将 会 对 物种 地 理 分 布 产 生 极 大 的 影响 。 

关于 物种 地 理 分 布 的 研究 方法 已 有 很 多 ,其 中 比较 热门 的 研究 方法 是 物种 分 布 模型 (Species Distribution 
Models,， SDMs) 。 物 种 分 布 模型 主要 是 利用 物种 的 分 布 数据 与 环境 数据 ,依据 特定 的 算法 估计 物种 的 生态 位 ， 
并 投影 到 景观 中 , 以 概率 的 形式 反映 物种 对 生境 的 偏好 程度 沾 。 在 所 有 物种 分 布 模型 中 ,最 大 粹 模型 
(Maxent) 相 比 而 言 预 测 效 果 最 好 ”| 。 

阔 叶 红 松林 树种 资源 丰富 , 红 松 (Pinus koraiensis ) 为 主要 成 林 树 种 ,并 与 其 他 针 痊 叶 树 种 如 鱼鳞 云 衫 
(Picea jezoensis) 、 奥 冷 杉 (Abies nephrolepis) 、 紫 要 (Tiia ammirensis)、 水 曲 柳 ( Frasinus mandshurica) 和 蒙古 标 
( Quercus mongolica ) 等 形成 针 阔 混交 林 '| ,是 东北 东部 山区 典型 的 地 带 性 植被 ,也 是 第 四 纪 冰 川 期 后 保存 完 
好 的 原始 裸子 植物 群落 ' ,其 分 布 中 心 在 东北 长 白山 、 张 广 才 岭 .小兴安岭 和 完 达 山 山 地 。 由 于 历史 上 的 过 
度 采 伐 , 阔 叶 红 松林 遭 到 严重 破坏 ,其 分 布 面积 急剧 萎缩 。 为 了 实现 阔 叶 红 松林 的 生态 恢复 及 可 持续 经 营 ,人 研 
究 阔 叶 红 松林 在 未 来 气候 变化 情况 下 的 潜在 分 布 区 显得 至 关 重 要 。 目 前 ,关于 阔 叶 红 松林 地 理 分 布 的 研究 还 
局 限于 单一 的 树种 ”中 ,而 单一 树种 的 分 布 区 域 不 能 直接 代替 阔 叶 红 松林 群落 的 分 布 区 域 ,并 且 这 些 研究 仅 
仪 探讨 了 为 数 不 多 的 环境 因子 ,如 温度 累积 指标 生长 度 : 日 (CDD) 和 水 分 指标 可 能 蒸 散 率 (PER ) 对 其 产生 
的 影响 。 本 文 以 阔 叶 红 松林 的 4 类 主要 建 群 树种 一 一 红 松 、 紫 要 水 曲 柳 和 蒙古 栎 为 研究 对 象 ,通过 最 大 
信 模 型 分 析 决 定 阔 叶 红 松 林地 理 分 布 的 主要 气候 和 地 形 因子 ,并 分 别 模拟 4 类 树种 的 潜在 分 布 区 ,根据 它们 
的 共同 分 布 区 域 来 决定 阔 叶 红 松林 的 分 布 区 ,最 后 利用 气候 变化 下 的 未 来 环境 数据 预测 阔 叶 红 松林 的 未 来 潜 
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在 分 布 区 ,为 阔 叶 红 松林 的 生态 恢复 和 经 营 管理 提供 理论 依据 。 
1 材料 与 方法 


1.1 研究 区 域 

阔 叶 红 松林 主要 分 布 于 我 国 东 北 地 区 ,其 地 处 亚 欧 大 陆 东 缘 , 地 理 坐 标 38"43' 一 53"23'N,118"50' 一 
135"05 中 ,包括 辽宁 吉林、 黑龙 江 3 省 以 及 内 蒙古 自治 区 东部 的 呼伦贝尔 市 .兴安 盟 .通辽 市 .赤峰 市 ,南北 长 
1600 多 km ,东西 宽 1400 多 km'"” , 南 临 黄海 和 渤海 ,东北 与 俄罗斯 接壤 ,东南 以 鸭绿江 为 界 与 朝鲜 民主 主义 
人 民 共 和 国 隔 江 相 望 ,西部 与 内 蒙古 自治 区 接壤 ' |。 
1.2 数据 采集 
1.2.1 树种 地 理 分 布 数据 

通过 中 国 数字 植物 标本 馆 ( http:/Awww.cvh.org.cn/ ) 教学 标本 资源 共享 平台 (http://mnh.scu.edu.cn/)、 
中 国 自然 保护 区 标本 资源 共享 平台 (http://www.papc.cn/ ) 等 数据 库 获得 标本 采集 地 信息 ,同时 查阅 《中 国 植 

志 》 和 东 三 省 植物 志 以 及 相关 研究 文献 确定 其 现在 分 布点 。 去 除 模糊 记录 的 分 布点 信息 ,对 具有 详细 信息 

的 分 布点 ,利用 百度 拾取 坐标 系统 (http://api.map.baidu.com/lbsapi/getpoint/index.html) 确 定 其 分 布点 的 经 纬 
度 。 另 外 ,为 提高 预测 的 可 靠 性 , 噜 除 引种 栽培 的 树种 分 布点 。 
1.2.2 气候 因子 

气候 基准 数据 来 源 于 WorldClim 数据 库 ( http://www.worldclim.org) ,该 数据 库 收集 了 1950 一 2000 年 全 球 
各 地 气象 站 的 每 月 气象 数据 ,采用 插值 法 生成 全 球 气 候 栅 格 数据 ,将 2000 年 作为 基准 年 。 本 文 基 于 ArcGIS 
10.0 平台 ,提取 了 东北 地 区 的 19 个 气候 因子 (Biol 一 Bio19) ( 表 1) ,空间 分 辨 率 为 30"( 约 1km?*)。 徐 影 所 
(2002) 对 国际 上 较 流行 的 5 种 大 气 环 流 模 式 (HadCM GFDL ECHAM CSIRO 以 及 CGCM) 进行 比较 研究 , 研 
究 结果 表明 ECHAM4 和 HadCM2 两 个 模式 对 东亚 和 中 国 地 区 的 气候 模拟 效果 最 好 ! ,本 文采 用 的 未 来 气候 
数据 由 大 气 环流 模型 HadCM 模拟 ,版 本 为 HadCM3。 研 究 涉及 的 3 种 排放 场景 来 自 《IPCC 排放 情景 特别 报 
告 》, 分 别 代 表 了 碳 排 放 高 速 增长 (SRES-A2) .中 速 增长 (SRES-A1B ) 和 低速 增长 (SRES-B1)。 其 中 A2 情景 
是 描述 了 一 个 非 均 衡 的 世界 :各 地 域 间 生产 力 方式 的 趋同 异常 缓慢 ,由 此 导致 人 口 持续 增长 ,经 济 发 展 主要 面 
向 区 域 , 人 均 经 济 增长 和 技术 变化 是 不 连续 的 ,并 低 于 其 他 情景 的 发 展 速 度 ;A1 情景 描述 的 世界 是 :经 济 增长 
迅速 全球 人 口 峰值 将 出 现在 本 世纪 中 叶 ,新 的 更 高 效 的 技术 被 迅速 引进 ,Al 情景 进一步 划分 为 3 组 情景 ,而 
AlB 是 能 源 使 用 较为 平衡 的 一 种 类 型 ;B1 情景 是 描述 了 一 个 趋同 的 世界 .全球 人 口 数量 与 Al 相同 ,但 经 济 
结构 向 服务 和 信息 经 济 方向 迅速 调整 , 伴 之 以 材料 密集 程度 的 下 降 ,以 及 清洁 和 资源 高 效 技术 的 引进 5 。 以 
上 3 种 排放 场景 3 个 未 来 时 段 (2020、2050、2080 年 ) 的 未 来 气候 数据 由 国际 热带 农业 中 心 ( CIAT, http:// 
www.ccafs-climate.org) 提供 ,空间 分 辨 率 为 30"(lkm”) 。 


表 1 研究 采用 的 环境 因子 


Table 1 Environmental variables used in the study 


类 型 代码 环境 因子 描述 
Type Code Environmental variables Description 
气候 因子 Climatic variable Biol” Annual mean temperature 年 平均 温度 
Bio2” Mean diurnal temperature range 每 月 最 高 温 与 最 低温 差 值 的 平均 值 
Bio3* Isothermality 等 温 性 
Bio4 Temperature seasonality 季节 性 温度 变异 
Bio5 Max temperature of warmest month 最 热 月 的 最 高 温度 
Bio6 Min temperature of coldest month 最 冷 月 的 最 低温 度 
Bio7 Temperature annual range 温度 的 年 较 差 
Bio8” Mean temperature of wettest quarter 最 湿 季 度 的 平均 温度 
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续 表 

类 型 代码 环境 因子 描述 

Type Code Environmental variables Description 
Bio9 Mean temperature of driest quarter 最 干 季 度 的 平均 温度 
Bio10 Mean temperature of warmest quarter 最 热 季 度 的 平均 温度 
Bioll Mean temperature of coldest quarter 最 冷 季 度 的 平均 温度 
Biol2”™ Annual precipitation 年 降雨 量 
Bio13 Precipitation of wettest month 最 湿 月 的 降雨 量 
Bio14 Precipitation of driest month 最 干 月 的 降雨 量 
Biol5™* Precipitation seasonality 降雨 季节 性 变化 
Biol6 Precipitation of wettest quarter 最 湿 季 度 的 降雨 量 
Bio17 Precipitation of driest quarter 最 干 季 度 的 降雨 量 
Biol8 Precipitation of warmest quarter 最 热 季 度 的 降雨 量 
Bio19 Precipitation of coldest quarter 最 冷 季度 的 降雨 量 

地 形 因子 Topographic variable ELE 上 Elevation 海拔 
SLO* Slope 坡度 
ASP” Aspect 坡 向 


* 经 筛选 后 输入 模型 的 因子 


12.3 地形 因子 和 底 图 

地 形 数据 来 源 于 中 国 科学 院 计算 机 网 络 信息 中 心 国 际 科学 数据 镜像 网 站 (http://www.gscloud.cn)。 下 
载 分 辩 率 为 90m 的 数字 高 程 地 图 ( DEM) , 从 中 提取 海拔 ,坡度 、 坡 向 ( 表 1) ,而 后 进行 重 采样 生成 分 辩 率 为 
30" 的 栅 格 数据 图 层 。 中 国 地 图 和 中 国 省 级 行政 区 划 图 来 自 国家 基础 地 理 信息 系统 网 站 (http://nfgis.nsdi. 
gov.cn/ )o 
1.3 数据 分 析 
1.3.1 树种 地 理 分 布 区 

将 树种 分 布点 整理 到 EXCFL 中 , 共 获 取 红 松 分 布点 167 个 ,蒙古 栎 分 布点 210 个 , 紫 根 分 布点 177 个 ,水 
曲 柳 分 布点 106 个 (图 1) 。 陈 新 美 等 对 Maxent 模型 预测 物种 分 布 所 需 样本 量 的 研究 表明 ,样本 量 大 小 对 
Maxent 模型 预测 物种 空间 分 布 的 精度 影响 不 大 "| ,所 以 本 文 数据 个 数 可 以 进行 模拟 。 
1.3.2 ”环境 因子 的 竹 选 

由 于 各 环境 因子 之 间 具 有 一 定 的 相关 性 ,如 果 直 接应 用 于 模型 ,可 能 会 产生 过 度 拟 合 现象 ,所 以 对 环境 因 
子 进 行 相关 性 分 析 后 才 可 以 应 用 于 物种 分 布 模型 71。 参 照 Yang 等 租 选 环境 因子 的 方法 上 ,对 22 个 环境 因 
子 图 层 进行 多 重 共 线性 分 析 (SPSS 19.0) 来 检验 图 层 之 间 的 相关 性 , 若 两 个 环境 因子 之 间 的 Pearson 指数 17| 
>0.8, 那 么 只 能 有 -一个 因子 选 和 模型 ,最 后 筛选 得 到 10 个 环境 因子 ,包括 7 个 气候 因子 ( Biol .Bio 2 .Bio 3、 
Bio 4、Bio 8、Bio 12 和 Bio 15) 和 3 个 地 形 因子 (海拔 ,坡度 . 坡 向 )( 表 1) 。 
1.3.3 物种 分 布 模型 

Maxent 模型 是 物种 分 布 模型 中 表现 较 好 的 一 种 模型 , 它 基于 可 以 免费 使 用 的 Maxent 软件 V3.3.3.k 
(http://www.cs.princeton.edu/ 一 schapire/maxent/) 。Maxent 软件 运行 输入 数据 包括 研究 区 域 的 一 组 环境 图 
层 和 一 个 物种 在 该 区 域 的 分 布 数据 。 参 照 Moreno 等 研究 方法 ' 中 ,随机 选取 25% 的 分 布点 作为 测试 集 (testing 
data) ,重复 运算 10 次 进行 建 模 , 其 他 参数 为 默认 值 。 采用 受 试 者 工作 特征 曲线 (receiver operating 
characteristic curve，ROC) 下 面积 ( Area under curve，AUC ) 来 评价 模型 的 拟 合 程度 。 选 择 AUC 值 最 大 的 一 次 
运算 结果 进行 物种 地 理 分 布 的 预测 , AUC 值 越 大 表明 拟 合 越 好 ,目前 以 AUC 值 0.5 一 0.6 为 较 差 ,0.6 一 0.7 为 
一 般 ,0.7 一 0.8 为 较 准确 ,0.8 一 0.9 为 很 准确 ,0.9 一 1 为 极 准确 。 在 环境 参数 设置 中 开启 刀 切 法 (Jackknife ) 来 
评价 各 个 环境 变量 的 权重 ,结合 每 个 环境 因子 的 贡献 百分比 来 判定 主导 因子 。 模 型 的 输出 格式 为 ASCLL 格 
式 文件 ,利用 ASCII to Raster 转换 为 栅 格 数据 , 按 栅 格 数值 的 大 小 将 研究 区 域 分 为 3 个 适宜 等 级 ,具体 参考 政 
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1 红 松 .蒙古 栎 、 紫 机 和 水 曲 柳 分 布点 的 地 理 位 置 


Fig.l Geographic locations of Pinus koraiensis, Quercus mongolica, Tilia amurensis and Fraxinus mandshurica 


府 间 气候 变化 专门 委员 会 (IPCC) 关 于 评估 可 能 性 的 划分 标准 !291 :P<0.05 为 不 适宜 区 ;0.05 大 P<0.33 为 低 度 
适宜 区 ;P 宇 0.33 为 高 度 适宜 区 。 本 文 将 4 类 树种 共同 具有 的 高 度 适宜 区 作为 阔 叶 红 松林 的 高 度 适宜 区 ,用 4 
类 树种 的 不 适宜 区 求 并 集 作为 阔 叶 红 松林 的 不 适宜 区 ,其余 作 为 低 度 适宜 区 。 


2 结果 与 分 析 


2.1 模型 模拟 的 准确 性 评价 

目前 ,ROC 曲线 分 析 法 已 经 广泛 应 用 于 准确 评价 物种 分 布 模型 的 预测 能 力 。 本 文 对 4 类 树种 的 分 布 进 
行 模拟 ,各 树种 的 AUC 值 如 下 ( 表 2)。 其 中 红 松 表现 最 好 ,其 次 是 水 遇 柳 、. 紫 机、 蒙古 标 。 当 Test AUC>0.8 
时 ,Maxent 模型 具有 很 好 的 预测 能 力 ,模拟 结果 保持 稳定 ,可 以 预测 物种 分 布 区 域 !'2 。 
2.2 影响 阔 叶 红 松 林 分 布 的 主导 环境 因子 

通过 对 4 类 树种 的 环境 因子 贡献 百分比 炙 加 分 析 ( 图 2) ,可 知 影响 阔 叶 红 松林 分 布 的 主导 环境 因子 是 年 
降雨 量 (Bio12) .季节 性 降雨 量 (Bio15) 海拔 (ELE) 、 年 平均 温度 ( Biol) 和 最 湿 季 度 的 平均 温度 ( Bio8) ,对 于 
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红 松 蒙古 栎 ,水 曲 柳 , 它 们 的 累计 贡献 率 都 大 于 85% ,对 于 紫 要 ,它们 的 累计 贡献 率 大 于 80% ,其 中 年 降雨 量 
为 最 主要 的 因子 。 基 于 Jackknife 模块 的 输出 结果 也 能 很 好 地 反映 这 些 因 子 对 各 树种 分 布 的 重要 性 。 
表 2 各 树种 模型 模拟 的 AUC 什 


Table 2 AUC values of Maxent models for given species 


树种 Species 曲线 下 面积 Area under curve 树种 Species 曲线 下 面积 Area under curve 
红 松 已 koraiensis 0.925 水 曲 柳 下. mandshurica 0.890 
紫 要 Tamurensis 0.859 蒙古 栎 Q. mongolica 0.847 


2.3 ”了 阀 叶 红 松林 适宜 分 布 区 的 空间 变化 趋势 分 析 
2.3.1 阔 叶 红 松林 潜在 分 布 区 的 模拟 

运行 Maxent 模型 后 输出 4 类 树种 的 潜在 分 布 区 ,利用 Reclassy 工具 划分 等 级 ,然后 辣 加 得 到 阔 叶 红 松林 
分 布 图 (图 3)。 阔 叶 红 松林 的 高 度 适宜 分 布 区 主要 分 布 在 长 白山 和 小 兴安 岭 地 区 , 东 起 吉林 田 春 市 黑龙江 
东 宁 县 , 北 界 位 于 黑龙 江 伊 春 市 . 稚 岗 市 更 北 县 ,向 南 经 小 兴安 岭 . 张 广 才 岭 .大 青山 .吉林 哈达 岭 .辽宁 本 溪 ， 
到 辽宁 册 岩 满族 自治 县 、 凤 城市 北部 、 桓 仁 满族 自治 县 。 阁 叶 红 松林 的 低 度 适宜 区 主要 分 布 在 东北 地 区 的 东 
南面 , 北 界 位 于 黑龙 江 黑 河 市 ,向 南 延伸 至 辽宁 普兰 店 市 ` 瓦 房 店 市 。 阔 叶 红 松林 的 高 度 适 宜 区 面积 占 研 究 区 
总 面积 的 11.69% , 低 度 适 宜 区 面积 占 研 究 区 总 面积 的 23% ,不 适宜 区 面积 占 研究 区 总 面积 的 65.31%( 表 3) 。 


SP Bio4 SLO ES Bio3 
Bio2 Biog ”BEEBBiol ELE 
Biol5 | Biol2 
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贡献 百分比 
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硬是 不 适宜 区 
[| 低 度 适宜 区 
[高 度 适 宜 区 


红 松 ” ” 紫 极 。 蒙古 村 水 曲 柳 0 


P. koraiensis T. amurensis OQ. mongolica F. mandshurica 


图 2 各 树种 的 环境 因子 贡献 百分比 图 3 阔 叶 红 松林 的 潜在 分 布 图 
Fig.2 Percentage contribution of environmental variables to the Fig.3 Potential distribution of broadleaved Pinus koraiensis forest 


species distribution models 


2.3.2 ”未 来 阔 叶 红 松林 潜在 分 布 区 的 预测 

由 基准 环境 数据 与 物种 分 布 数据 建立 的 模型 ,预测 效果 良好 。 在 此 基础 上 ,将 物种 和 环境 因子 之 间 的 联 
系 投影 到 未 来 气候 情景 中 用 于 预测 4 类 树种 的 未 来 分 布 ,得 到 阔 叶 红 松林 在 3 种 气候 情景 下 ,2020、2050、 
2080 年 的 预测 分 布 图 (图 4) 。 

在 A2、AlB .Bl1 气候 情景 下 , 阔 叶 红 松林 高 度 适宜 区 将 主要 集中 分 布 于 长 白山 自然 保护 区 和 小 兴安 岭 部 
分 地 区 。 其 分 布 南 界 与 北 界 都 向 北 移动 ,在 A2、A1B 情景 下 的 2080 年 尤为 明显 。 

在 A2 气候 情景 下 , 阔 叶 红 松林 的 高 度 适 宜 区 面积 波动 不 定 ,由 原来 的 11.69% 怠 减 到 2020 年 的 4.77%， 
到 2050 年 增加 到 5.66% ,而 后 2080 年 又 减少 为 0.60% ; 益 叶 红 松 林 的 低 度 适宜 区 面积 从 2000 年 的 23 多 呈现 
增加 的 趋势 ,到 2080 年 为 37.18% (图 4, 表 3)。 
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图 4 不 同 气候 情景 下 , 阔 叶 红 松林 的 潜在 分 布 预测 图 


Fig.4 Predicted potential distribution of broadleaved Pinus koraiensis forest under different climate change scenarios 


1 : A2 气候 情景 ; 工 : A1B 气候 情景 ; 于: B1 气候 情景 ;1: 21 世纪 20 年 代 ;2: 21 世纪 50 性 


FE 代 ;3: 21 ] 


情景 下 21 世纪 50 年 代 


此 纪 80 年 代 。 如 : I 一 2, 指 AlB 气候 


在 AlB 气候 情景 下 , 阔 叶 红 松林 的 高 度 适宜 区 面积 表现 为 减少 趋势 ,由 基准 年 的 11.69% 又 减 到 2020 年 
的 4.71% ,到 2050 年 减 小 到 2.17% , 而 2080 年 缩 少 至 0.46% ; 阔 叶 红 松林 的 低 度 适宜 区 面积 从 2000 年 的 23% 


增加 到 2020 年 的 28.01% ,至 此 接近 稳定 。 
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表 3 阔 叶 红 松林 的 不 同 适 宜 等 级 面积 百分比 
Table 3 Habitat suitability classes area percent of broadleaved Pinus koraiensis forest 


面积 百分比 例 Percentage of area/ % 


气候 变化 情景 年 份 


Climate change scenario Year 不 适宜 区 低 度 适宜 区 高 度 适宜 区 
Unsuitable area Low suitable area High suitable area 

基准 current 2000 65.31 23.00 11.69 

A2 气候 情景 A2 scenario 2020 70.88 24.35 4.77 
2050 61.91 32.44 5.66 
2080 62.22 37.18 0.60 

A1B 气候 情景 A1B scenario 2020 67.27 28.01 4.71 
2050 69.26 28.57 2.17 
2080 71.32 28.22 0.46 

B1 气候 情景 B1 scenario 2020 63.32 23.54 13.14 
2050 63.09 30.10 6.81 
2080 65.30 30.34 4.36 


而 在 Bl 气候 情景 下 , 阔 叶 红 松林 的 高 度 适宜 区 面积 开始 稍 有 增加 而 后 减少 ,由 11.69% 增 长 至 2020 年 的 
13.14% ,但 在 2050 和 2080 年 其 面积 分 别 降 为 6.81% 和 4.36%; 阔 叶 红 松林 的 低 度 适宜 区 面积 从 2000 年 
(23% ) 到 2020 年 (23.54% ) 几乎 没有 变化 ,但 是 到 2050 年 增加 为 30.10% , 而 后 保持 稳定 。 


3 讨论 


已 有 研究 表明 ,限制 植物 地 理 分 布 的 主要 因子 是 植物 的 耐寒 性 .完成 生活 史 所 需 的 热量 供应 和 可 利用 
水 加 。 笑 晓 洁 等 研究 表明 影响 蒙古 栎 的 最 主要 因子 是 年 均 降 水 量 :“ ,本 文 得 到 与 其 相同 的 结论 。 孙 晓 红 研 
究 表明 影响 红 松 地 理 分 布 的 最 主要 气候 因子 是 最 冷 月 最 热 月 平均 温度 差 ,其 次 是 年 降水 量 !” ,与 本 研究 结果 
有 一 定 的 差异 。 其 主要 原因 是 不 同 的 研究 所 采用 的 环境 因子 不 同 ,本 文 增加 了 地 形 因子 作为 环境 因子 ,此 外 ， 
由 于 环境 因子 的 筛选 过 程 不 同 ,也 会 将 不 同 的 环境 因子 选 人 模型 。 

对 东北 地 区 气温 和 降水 变化 的 研究 表明 ,在 A2、A1B .BI 三 种 排放 情景 下 未 来 100 年 的 气温 和 降水 总 体 
趋势 均 呈 逐渐 增加 的 趋势 。 此 外 ,在 A2 情景 下 气温 和 降水 的 增幅 最 高 ,其 次 是 AlB 情景 ,变化 幅度 最 小 的 是 
B1 排放 情景 .1 。 有 研究 表明 红 松 地 理 分 布 区 面积 与 降水 量 呈 正 相关 关系 ,与 年 平均 温度 呈 负 相关 关系 '*| ， 
而 全 球 气候 变化 将 导致 东北 地 区 温度 明显 升 高 ,温度 带 北 移 , 降 水 量 有 所 增加 ,但 抵消 不 了 由 温度 升 高 造成 的 
蒸 散 量 的 增加 ,因此 全 球 气 候 变 化 将 使 东北 东部 山区 面临 暖 干 的 气候 变化 趋势 ,所 以 实际 上 气候 的 暖 干 化 促 
使 阔 叶 红 松林 分 布 区 面积 有 减少 的 趋势 ,生态 适宜 性 显著 下 降 情 ; 。 程 肖 侠 等 研究 表明 ,在 气候 变 暖 降水 增加 
的 背景 下 , 红 松 有 成 为 大 兴安 岭 地 区 森林 主要 树种 的 趋势 :2 ,而 本 文中 在 3 种 气候 情景 下 , 低 度 适宜 区 都 有 
显著 的 北 移 ,到 2020 年 ,其 分 布 区 域 扩展 到 大 兴安 岭 地 区 ,在 A2、A1B 气候 情景 下 到 2080 年 ,大 兴安 岭 的 部 
分 地 区 也 将 会 成 为 阔 叶 红 松林 的 高 度 适 宜 区 。 但 对 于 小 兴安 岭 地 区 的 研究 却 有 所 不 同 , 当 未 来 年 降水 增加 
10% 左 右 \ 年 气温 增幅 大 于 5%C 时 , 林 徐 模型 的 模拟 结果 表明 阔 叶 红 松林 将 被 蒙古 栎 、 紫 要 和 和 裂 叶 榆 为 主 组 成 
的 阔 叶 林 所 取代 "3 ,在 CGCM2 情景 下 ,也 能 够 得 到 相同 的 结论 55 ,而 本 文选 取 SRES-A2、SRES-A1B、SRES- 
B1 三 种 气候 情景 ,结果 有 很 大 的 差别 。 原 因 是 物种 分 布 模型 目前 有 很 多 不 确定 性 ,首先 ,在 模型 选取 上 ,现在 
有 很 多 种 类 的 物种 分 布 模型 ,不 同 的 物种 分 布 模型 预测 结果 有 所 差异 ,其 次 ,环境 因子 的 选取 也 有 很 大 变化 ， 
选取 不 同 的 气候 情景 ,预测 结果 会 有 所 不 同 '* 。 此 外 ,在 A2 气候 情景 下 ,从 2020 年 到 2050 年 阔 叶 红 松林 的 
高 度 适宜 区 面积 增加 ,其 原因 是 阔 叶 红 松林 的 适宜 性 等 级 分 布 是 采取 4 类 树种 的 到 加 方式 得 到 的 ,因此 可 能 
会 造成 累计 误差 。 

物种 分 布 区 是 物种 生态 与 进化 历史 的 复杂 表达 形式 ,在 不 同 的 时 空 尺 度 上 受到 多 种 因素 ,不同 强度 的 控 
制 .” 3 。 通 常 认 为 决定 物种 分 布 区 共有 四 种 因素 :环境 因素 .生物 因素 物种 的 扩散 能 力 .物种 适应 新 环境 的 
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进化 能 力 ' ”1 。 而 本 文 仅仅 考虑 环境 因素 中 的 部 分 环境 因子 作为 预测 指标 进行 预测 ,与 实际 还 有 一 定 的 差 
距 , 还 需要 做 进一步 的 研究 。 


目前 , 阔 叶 红 松林 主要 分 布 于 长 白山 和 小 兴安 岭 地 区 。 但 是 在 A2、A1B 、B1 气候 情景 下 预测 阔 叶 红 松林 
的 分 布 ,结果 有 很 大 的 差异 ,其 分 布 界线 北 移 、 高 度 适宜 分 布 区 域 面 积 减 小 ,其 中 长 白山 地 区 的 阔 叶 红 松林 高 
度 适 宜 分 布 区 更 集中 于 长 白山 自然 保护 区 。 此 外 ,气候 变化 程度 越 严 重 ,对 阔 叶 红 松林 的 分 布 影响 也 越 大 , 尤 
其 是 在 A2、A1B 气候 情景 下 ,到 2080 年 阔 叶 红 松林 的 高 度 适宜 区 面积 所 占 研究 区 域 面 积 不 足 1% ,说 明 如 果 
不 加 以 保护 环境 而 使 气候 得 以 改善 ,那么 阔 叶 红 松林 在 未 来 有 可 能 逐渐 消失 于 东北 地 区 。 而 影响 阔 叶 红 松林 
分 布 的 主导 环境 因子 主要 是 年 降雨 量 (Bio12) .季节 性 降雨 量 (Bio15) 海拔 (ELE) .年 平均 温度 (Biol) 和 最 湿 
季度 的 平均 温度 (Bio8 ) 。 
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